Određivanje uopštene kompleksne funkcije
assume(0 < R & 0 < C & 0 < L )
zamena = L == C*R^2
zamena = 
vrednosti = [C == 1, R == 1]
vrednosti = 
Z1 = R + s*L
Z1 = 
Z2 = 1/(s*C + 1/R)
Z2 =

H(s) = subs(Z2/(Z1 + Z2), lhs(zamena), rhs(zamena))
H(s) =

simplify(H(s))
ans =

Amplitudska i fazna karakteristika
Hjw(w) = simplify(subs(H(s), s, 1i*w))
Hjw(w) =

Aw(w) = simplify(abs(Hjw(w)))
Aw(w) =

phi(w) = expand(angle(Hjw(w)))
phi(w) =

Crtanje amplitudske i fazne karakteristike
fplot(w, subs(Aw(w), lhs(vrednosti), rhs(vrednosti)), [0,4])
fplot(w, subs(phi(w), lhs(vrednosti), rhs(vrednosti)), [0,4])
Određivanje propusnog opsega
Aref = Aw(0)
Aref =

w3dB = solve(Aw(w) == Aref/sqrt(2), w)
w3dB =

BandPass3dB = [0, w3dB]
BandPass3dB =

wg1 = subs(w3dB(1), lhs(vrednosti), rhs(vrednosti))
wg1 = 
Označavanje propusnog opsega u grafiku
fplot(w, subs(Aw(w), lhs(vrednosti), rhs(vrednosti)), [0,4])
rectangle('Position',[0, 0, double(wg1), double(Aref/sqrt(2))],...
'FaceColor',[1 0.9 0.9],...